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P ; ABSTRACT 



Using a thermal-chemical model for the generic T-Tauri disk of D'Alessio 
et al. (1999), we estimate the strength of the fine-structure emission lines of Nell 
and Nelll at 12.81 and 15.55 /im that arise from the warm atmosphere of the 
disk exposed to hard stellar X-rays. The Ne ions are produced by the absorption 
of keV X-rays from the K shell of neutral Ne, followed by the Auger ejection of 
several additional electrons. The recombination cascade of the Ne ions is slow 
because of weak charge transfer with atomic hydrogen in the case of Ne"'"^ and 
by essentially no charge transfer for Ne"*". For a distance of 140 pc, the 12.81 /im 
line of Nell has a flux ~ 10~^'^ergcm~^s~^, which should be observable with 
the Spitzer Infrared Spectrometer and suitable ground based instrumentation. 
The detection of these fine-structure lines would clearly demonstrate the effects 
of X-rays on the physical and chemical properties of the disks of young stellar 
objects and provide a diagnostic of the warm gas in protoplanetary disk atmo- 
spheres. They would complement the observed H2 and CO emission by probing 
vertical heights above the molecular transition layer and larger radial distances 
that include the location of terrestrial and giant planets. 
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Subject headings: accretion, accretion disks — infrared: stars — planetary sys- 
tems: protoplanetary disks — stars: formation — stars: pre-main sequence — 
X-rays: stars 

1. Introduction 

Young stellar objects (YSOs) are observed to be strong emitters of X-rays (e.g., Feigelson 
& Montmerle 1999). The X-rays have the potential to affect the physical properties of nearby 
circumstellar material (e.g., Glassgold, Feigelson & Montmerle 2000) and thus the course of 
star formation (Silk & Norman 1980; Pudritz & Silk 1987). Tsujimoto et al. (2005) observed 
X-ray fluorescence from cold Fe ions in seven sources in the Orion Nebula Cluster, which they 
have interpreted to be the result of absorption of hard stellar X-rays by cool circumstellar 
disk material. This confirmation of the interaction of stellar X-rays with the disks of YSOs 
raises the question of whether there are specific diagnostics for the effect of X-rays on the 
physical and chemical properties of these disks, which are believed to play a crucial role in 
the formation of stars and planets. Although the gas is the main reservoir of mass for much 
of the lifetime of a YSO disk, much less is known about the gas than the easier to measure 
dust. Thus we address the issue of X-ray effects in the broader context of the search for 
gaseous diagnostics of disk gas. 

This type of question has a long history in studies of the interstellar medium of our own 
and external galaxies, where definitive demonstrations of the effects of X-rays have been 
elusive because other external heating and ionizing radiations, such as UV radiation and 
cosmic rays, can produce similar effects. Some of the diagnostic possibilities discussed in 
the literature are multiply charged ions (Dalgarno 1976; Langer 1978); specific molecules, 
especially ions and radicals (Krolik & Kallmann 1983; Neufeld, Maloney & Conger 1994; 
Lepp & Dalgarno 1996; Sternberg, Yan & Dalgarno 1997; Yan & Dalgarno 1997); near- 
infrared H2 ro-vibrational transitions (Lepp & McCray 1983; Draine & Woods 1989; Gredel 
& Dalgarno 1995; Tine et al. 1997); and fine-structure fines (Maloney, HoUenbach & Tielens 
1996). 

Many of these possibilities are now being re-investigated in the context of chemical 
modefing of the disks and envelopes of YSOs irradiated by X-rays and UV radiation (e.g., 
Aikawa & Herbst 1999; Markwick et al. 2002; Gorti & HoUenbach 2004; Millar & Nomura 

2005; Stauber et aL2005). An important aspect of the X-ray irradiation of protoplanetary 
disks is that the upper atmosphere close to the YSO is heated to relatively high tempera- 
tures ~ 5,000K (Glassgold, Najita & Igea 2004, henceforth GN104; Alexander, Clarke & 
Pringle 2004). The temperature then decreases going down towards the disk midplane. The 
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accompanying thermal-chemical structure consists of a hot partially-ionized atomic layer on 
top of a cool molecular layer with a much smaller ionization fraction. In between, there 
is a warm partially-molecular layer with gas temperatures ~ 500 — 2, 000 K. This layered 
structure of X-ray irradiated disks can be traced back to earlier studies of various astrophys- 
ical environments, e.g., Halpern & Grindlay (1983), Krolik & Kallman (1983) and Lepp & 
McCray (1983). The thermal properties of the disk atmosphere above the cool molecular 
layer play a critical role in determining the observational signatures of disks. 

The potential for ambiguity in the identification of X-ray diagnostics also arises in 
protoplanetary disks because of the presence of stellar UV radiation which can heat the 
gas in the surface layers via the photoelectric effect on small grains and polycyclic aromatic 
hydrocarbons compounds or PAHs (Kamp & DuUemond 2004; Nomura & Millar 2005). In 
the far outer regions of disks, interstellar UV radiation and galactic cosmic rays may also 
play a role in heating and ionizing the gas. We suggest that the mid-infrared fine-structure 
lines of Nell at 12.81 /xm and of Nelll at 15.6 /im may provide unique signatures of X-ray 
ionization of the disks of low-mass YSOs. 

Neon fine-structure lines have proved invaluable for studying ionized regions in our own 
as well as external galaxies (e.g., see the recent observations by Jaffe et al. 2003 of ultra 
compact HII regions in Monoceros R2, by Sturm et a/. 2001 of AGN, by Dcvost et a/. 2004 
of the starburst galaxy NGG253 and by Weedman et a/. 2005 of Scyfcrt galaxies). The 
Ne ions in HII regions and galaxies are generated by the Lyman continuum photons from 
massive stars and by X-rays. Generating these ions with UV photons is much more difficult 
in the environment of low-mass young stellar objects. Cosmic rays are also unlikely to be 
effective for this purpose because they are excluded by the strong winds generated by YSOs 
(Glassgold, Najita & Igea 1997). 

Low-mass YSOs arc known to be strong emitters of moderately hard keV X-rays (e.g., 
Feigelson & Montmerle 1999), and we suggest that the Ne ions are mainly produced by X- 
rays. The basic idea is that the high ionization potentials of Ne and Ne+ (21.56 and 41.0 eV, 
respectively) preclude the production of these ions by the FUV radiation (A > 911.6 A) 
emitted by T Tauri stars. We will show that the hard X-rays emitted by YSOs are effective 
in generating Nc+ and Nc^+ over a range of radii in the atmospheres of the disks of these 
stars. The emission arises in the same region of the disk from which photo-evaporation is 
believed to originate. Thus the Ne ion fine-structure line emission can serve as a diagnostic 
of the source of an evaporative fiow as well as the signature of X-ray irradiation. 

T Tauri stars are also likely to emit significant amounts of EUV radiation (A < 911.6 A ) 
(e.g., Alexander, Clarke and Pringle 2005), although it is much less well characterized than 
their X-ray emission. EUV radiation has been invoked to photo-evaporate disks (reviewed by 
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HoUenbach et al. 2000, DuUemond et al. 2006). However, EUV photons are rapidly absorbed 
by atomic hydrogen in the partially- ionized wind and disk of YSOs (Alexander, Clarke and 
Pringle 2003) so that, except for a fully- ionized region near the star (Hollenbach et al. 1994), 
the EUV flux that reaches substantial radii and depths in disks to generate Ne ions is 
probably very small. The small ionized region generated by stellar EUV radiation may be 
the source of some Ne fine-structure line emission (U. Gorti and D. Hollenbach, private 
communication 2006). 

The plan of this article is as follows. First, we lay the foundation for calculating the 
fluxes of the X-ray generated Ne flnc-structurc emission by developing the ionization theory 
for the Ne ions. We take into account production by X-rays and destruction by charge 
exchange as well as radiative recombination in a form appropriate to the thermal-chemical 
structure of X-ray irradiated protoplanetary disk atmospheres. In §3, we calculate the fluxes 
of the fine-structure lines for the disk atmosphere of a typical T Tauri disk, where the electron 
densities are generally sub-critical. In §4 we discuss the observability of these lines, their 
potential astrophysical significance, and the uncertainties in the theory. 



2. Ionization Theory 

We start by discussing the production and destruction of Ne ions in the inner region of a 
protoplanetary disk atmosphere, within tens of AU. We consider only X-rays for ionizing the 

ncTitral and low-ionization forms of the Ne atom, since FUV photons have insufficient energy 
to ionize them, EUV is absorbed over short distances in mainly neutral matter, and cosmic 
rays are largely excluded from the inner regions, as mentioned in §1. Collisional ionization 
by thermal electrons is unimportant for the relevant temperatures, T < 10, 000 K. 

The result of an X-ray absorption by a low-ionized Ne ion depends on whether the X-ray 
energy E is less than or greater than the K-shcll threshold, = 0.870, 0.903, . . . keV for 
Nel, Nell ... . According to simple one-electron theory (e.g., Kaastra & Mcwc 1993), 
an L-shell electron is ejected when E < E^, increasing the ion charge by one. For E > 
Ek, the cross section for L-shell photoionization is small compared with the production 
of a K-shell vacancy, and the latter leads mainly to the ejection of another electron by 
the Auger process. According to one-electron theory (ignoring correlations), the branching 
ratio for single ionization is only 2% (Kaastra & Mewe 1993). The actual situation is more 
complicated, especially with regard to the production of more than two electrons, as was 
shown in laboratory experiments on neutral Ne by Carlson & Krause and collaborators 
(Krause et al. 1964; Carlson & Krause 1965a,b). First of all, there is a small but significant 
probability, ~ 10%, that L-shell photoionization ejects two electrons. Second, the probability 
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that K-shell vacancies produce just one electron is somewhat larger (~ 6%) than predicted by 
one-electron theory ignoring electron correlations. Third, and most important, ~ 30% of K- 
shell vacancies produced in the photoionization of neutral Ne lead with significant probability 
to higher ions such as Ne^"*" and Ne*^+. However, in partially-ionized regions characteristic 
of the atmospheres of protoplanetary disks , where atomic hydrogen is the dominant species 
in contrast to fully ionized H II regions, Ne^"*" is rapidly destroyed by charge exchange with 
atomic hydrogen, 

Ne^+ + H^Ne^+ + H+. (2-1) 

The cross section for this reaction has been measured to be large and rising down to 0.1 eV by 
Rejoub et al. (2004); it corresponds to rate coefficients for charge transfer ~ 10~^cm^s~^. In 
essence, ions beyond Ne^+ that arise from K-shell vacancies in Ne atoms and low-ionization 
ions are rapidly converted back to Ne^"^ by one or more charge transfers with atomic hydro- 
gen. 

Charge exchange with atomic hydrogen is also important in the ionization balance of 
Ne, Ne+ and Ne^+. Charge transfer of Ne^+ to H was studied by Dalgarno and collaborators 
25 years ago (Butler, Bender & Dalgarno 1979; Dalgarno, Butler & Heil 1980; Butler, Heil & 
Dalgarno 1980; Butler & Dalgarno 1980). Their conclusion was that the rate coefficient for 
this process is exceedingly small, but that radiative charge exchange takes place with rate 
coefficient k2 ~ lO^^^cm'^s^^. Although this is much smaller than the (essentially Langevin) 
rate coefficient for fast charge transfer ('^ 10~^cm^s~^), it cannot be ignored in partially- 
ionized regions because its effect is the same order of magnitude as radiative recombination, 
oc Xea2, where 0:2 is the recombination rate coefficient for Ne^"*" and is the fractional 
ionization. There are no reports on charge transfer from Ne"*" to H, but the potential en- 
ergy curves calculated for this system by Pendergast, Heck and Hayes (1994) indicate that 
radiative charge exchange is highly forbidden (A. Dalgarno, private communication). We 
conclude that Ne"*" is primarily destroyed by radiative recombination and, to a lesser extent, 
by X-ray photoionization into higher ions. 

These considerations on Ne+ and Ne^+ charge transfer to H provide the basis for a simple 
set of steady-state ionization balance equations. The steady state assumption is justified by 
the high electron density in the region relevant for the emission of the fine-structure lines, 
typically ~ lO^cm"^. At these high densities, the time scale for the slowest process, radiative 
recombination, is ~ 1 yr, the same order or less than the orbital time and, more important, 
significantly less than the accretion time scale. We denote the fraction of Ne in charge state 
n by (n = 0, 1, 2). We express the difference between X-ray ionization below and above 
the Nel K-edge by writing the total X-ray ionization rate as the sum of a soft and a hard 
component, 

C(n) = CH + aW (n = 0,l). (2-2) 
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FoUowing the above discussion, we introduce branching ratios Bi and B2 for soft X-ray 
ionization of Nc to produce Ne"*" (Bi 0.9) and for hard X-ray ionization of Nc to produce 
Ne^"*" {B2 ~ 0.94, when fast charge exchange of higher ions back to Ne^"*" is included). The 
steady balance equations for Ne+ and Ne^"*" are then: 



Bi 



00) 



(1-52) 



0(0) 



Xo + {a2Xe + k2)x2 



01) 



+ {aiXe + ki) 



5oao) + ^,c.(o) 



Xo H Xi = (a2Xe + k2) X2, 



(2-3) 



(2-4) 



where ub. is the number density of hydrogen nuclei. These equations can easily be solved in 
terms of the ion abundance ratios. 



Xi 
Xo 



X2_ 
Xo 



B, 



a(o)/nH 



1-^1 



CKiXe + k\ ' 

C(0)/nH 



C(l)/nH Xi 



(2-5) 



(2-6) 



■ 0:2X0 + k2 a2Xe + k2 a2Xe + k2 Xo 

with the proviso that Xo + Xi + X2 = x^^, where is the disk abundance of gaseous Ne. 
When we recall from §2 that ki is neghgible and ignore the small quantities (1 — S„) and 
the small difference between ({0) and C(l)) we obtain the following relations: 



Xj_ 
Xo 



C(Ne)/nH 



X2 

Xo 



aiXe 

a(Ne)/nH , C(Ne)/nH Xi 



(2-7) 



(2-8) 



a2Xc + k2 a2Xc + k2 xq 
where C(Ne) = (^s(Ne)-|-C/i(Ne) is the ionization rate for neutral Ne per Ne nucleus. The direct 
X-ray ionization rates for Ne can be calculated from the absorption cross section (Morrison 
& McCammon 1983; Wilms et al. 2000). Here we adopt an approximate treatment in which 
the Ne rates are scaled from the ionization rate per H nucleus, according to the ratio 
of the Ne cross section to the mean cross section that defines (. Following GNI04, we 
assume that the X-ray luminosity and spectral temperature are Lboi = 2 x 10^°ergs~^ and 
kTx = 1 keV, so that the X-ray spectrum extends from several hundred eV to a few KeV. 
Good representative figures before attenuation are C(Nc) ~ 40^ and (^;j(Ne)/(^(Ne) ^ 7/8. 
These ratios are actually a function of position because the X-ray spectrum hardens due 
to the preferential absorption of soft X-rays. Thus the ratio C/i(Ne)/^(Ne) tends to unity 
at large column density. The ratio C(Ne)/C is much greater than one because the X-rays 
emitted by young stellar objects have characteristic energies of the same order or larger 
than the Ne K-edge near 0.9 keV, where the Ne absorption cross section is maximum. In 
particular, the photoionization cross section of Ne for keV X-rays is much larger than the 
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cross section summed over a cosmic abundance mix (Morrison and McCammon 1983; Wilms 
et al. 2000) that determines the ionization rate per hydrogen nucleus. If the stellar X- 
ray spectrum contained few keV photons, the ratios C(Ne)/C and Ch(Ne)/C(Ne) would be 
very much smaller than the representative values quoted above. The secondary electrons 
generated by the primary photoionization can ionize hydrogen as well as a heavy element 
like Ne. But the electronic ionization cross section for Ne is generally no more than twice 
those for hydrogen and helium. Following the discussion of Maloney, Hollenbach & Tielens 
(1996), we estimate that the contribution of the secondary electrons to the X-ray ionization 
rate of Ne is < 2^, and thus small compared to the direct X-ray ionization. 

We will find in §3 that, for the regions of the disk responsible for the emission of the 
fine-structure lines, the fractions Xi/xq and X2IXQ are usually < 0.1 and thus 

xi _ C(Ne)/nH _ X2 ^ Cfe(Ne)/nH _ . . 

— = Ol, r ; ■ = 02. l^-yj 

Xq tti^e Xq ^2 + Q^2 Xq 

Because 02 > ai and /c2 is finite though small, these equations tell us that there is substan- 
tially more Ne+ than Ne2+. 

Over most of the region of interest, the electron fraction is approximately given in terms 
of the total ionization rate by the simple formula. 



c 



nH«(H+) 



1/2 

(2-10) 



where q;(H+) is the total recombination rate of H+. Therefore, according to equations 12-91 
and I2-10[ the fractions of Ne in the first two ion states depend on different powers of the 
ionization parameter C/^h and, equivalently, the electron fraction, 

- ~ iC/nnf" ~ Xe, ^ ~ C/nn ~ xl (2-11) 

Xq Xq 

The actual ionization fractions used in the calculations of §3 are, 

/I ai {ai + h)a2 , . 

xo/x^e = 2^ xi/xNe = a;2/a;Ne = 2 ' (2-12) 

where = 1 + Oi + (oi + 6)03 with h = C/i(Ne)/C(Ne) ~ 0.875. 

In addition to their utility in the numerical calculations, these analytic formulae provide 
some insight into the physics of Neon ionization in X-ray irradiated protoplanetary disk 
atmospheres. In the first member of Eq. 12-lH the approximate proportionality of xi/xq to 
the electron fraction is similar in form to the well known ratio of 0+ to O when O"*" and 

are strongly coupled by fast, near-resonant, charge exchange (e.g., Osterbrock 1987). In 
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the present case, Ne"*" and H"*" are not connected directly by charge exchange, but they are 
hnked by both being produced by X-rays. A semi-quantitative expression of this relation, 
based on Eq. 12-11^ is, 

^(Ne+) a(H+)C(Ne) 

Because the ratio C(Ne)/C ^ 1 and ai is not very much smaller than q;(H"'"), a;(Ne^)/a;(Ne) 
is much larger than x^- Although Eq. 12-131 expresses much of the essential physics of the 
X-ray ionization of Ne for protoplanetary disks, it is inadequate for quantitative estimates 
because it breaks down for very small and very large Xe- 

In the next section, we use the detailed Ne ionization theory presented above to estimate 
the Ne ion line fluxes for a T-Tauri disk. We adopt, for this purpose, the thermal-chemical 
disk structure of an existing model that includes stellar X-ray heating and ionization (GNI04) 
but which did not include the Ne ionization theory discussed above. The GNI04 model allows 
for the inclusion of additional mechanical heating of the disk atmosphere, e.g., from accretion 
or wind-disk interaction, using the following formula for the volumetric heating rate 

9 

Tmech = -ahpc^VL, (2-14) 

where an is a phenomenological parameter. Figure 1 is an example of the results of GNI04 
at a radius of i? = 5 AU; additional results for a range of radii are presented below in Figures 
2 and 3. Figure 1 compares the vertical variation of the gas and dust temperatures. The 
abscissa represents the height above the midplane z in terms of the vertical column density 
A'h of hydrogen from oo down to z. The calculations are based on the D'Alessio et al. (1999) 
fiducial dust model for a T Tauri star with mass M^, = 0.5Mq, radius = 2Rq, and effective 
temperature T* = 4000 K. The accretion rate is M = 1O~^M0 yr~^, and the disk thickness 
at 1 AU is ~ 100 g cm~^ and varies elsewhere roughly as 1/R. Heating of small grains by 
the stellar optical and near infrared radiation produces an inversion in the dust temperature 
Td over the midplane value, as shown by the dashed curve. X-ray heating generates a much 
larger inversion in the gas temperature Tg (solid line labeled ah = 0.01) in the region above 
the warm dust. When mechanical heating is included (the solid curve labeled at = 1-0), the 
hot gas extends even further down into the disk. The gas and dust temperatures approach 
one another near the midplane (large iVn), where the volume densities tt-h are big enough to 
ensure good thermal coupling between the gas and the dust. In between the hot gas layer at 
the top of the atmosphere and the relatively cool midplane region is a warm transition zone 
(iVn > lO^^cm"^) where Tg ?a 500 — 2000 K. Gas temperature inversions can also be obtained 
with UV heating of small grains and PAHs (Kamp & Dullemond 2005, Nomura & Millar 
2005). Major chemical transitions occur between 10^^ and lO^^cm"^. Significant levels of 
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H2 and full formation of CO occur near A^h ~ 1 — 2 x lO^^cm ^, and complete formation of 
H2 is achieved at A^h = 8.5 x 10^^ cm~^. 

Figure 2 shows gas temperature profiles for radii ranging from R — 1 — AO AU for the 
case afi — 0.01. They exhibit a hot gas layer on top of a cool one, with temperatures 
generally decreasing with increasing radial distance. The temperature at the top of the 
atmosphere decreases with radius due to inverse square dilution of the stellar X-rays. The 
temperature near the midplane decreases with distance because the viscous heating in the 
D'Alessio model decreases. The results for the case where mechanical heating dominates are 
similar, except that the transition from hot to cold occurs deeper down in the atmosphere. 
Beyond R = 20 AU, there is a qualitative change in the temperature profiles in that the 
maximum temperature drops from several thousand to several hundred K. Figure 3 shows 
clearly how X-ray irradiation can warm disk atmospheres out to large radii of the order of 
tens of AU. It is also noteworthy that the near-midplane temperature at large radii can fall 
below the grain freeze-out temperature for many volatile species. Although the freeze-out 
process is omitted from our present chemical model, it does not affect the calculation of the 
Ne fine-structure emission, which arises from the upper atmosphere of the disk where the 
temperature is high and freeze-out does not occur. Regions where freeze-out can occur do 
not contribute to the emission because the upper levels of the transitions, with excitation 
temperatures ~ 1000 K, are difficult to excite with electron coUisions at low temperatures. 

Figure 3 shows that the electron fraction decreases smoothly at high altitudes. It then 
drops more rapidly to a level ~ 10^^—10^^, at a depth close to where Ty begins to approach 
Td and near where carbon becomes almost fully associated into CO. Figure 3 demonstrates 
that stellar X-rays can generate substantial levels of ionization in disk atmospheres out 
to large radii. The variation in gas temperature and ionization in Figures 2 and 3 are 
directly relevant to the strength of the Ne line fluxes. Above the thermal-chemical transition, 
especially in the range near A^h = 10^^ — 10^°cm~^, the warm ionized conditions are conducive 
to obtaining large Ne ion abundances and significant populations of the upper levels of the 
Nell 12.81 /xm and Nelll 15.55 /xm transitions. 

3. Flux Calculation 

We estimate the emission from the Ne I and Ne II fine-structure lines from a face-on disk 
in the optically-thin approximation. The lines originate in magnetic dipole transitions from 
the Nell doublet and the Nelll triplet, whose properties are given in Table 1. 
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Table 1: Ne Fine Structure Levels 



Ion 


A (/xm) 


Ju — Jl 


Tui{K) 


Aui (s-i) 


T-^/^ncr{ul) (cm-3) 


Auihuui (erg s ^) 


Ne+ 


12.81 


1/2-3/2 


1122.8 


8.59 X 10-3 


5.53 X 103 


1.332 X 10-^5 


Ne2+ 


15.55 


1-2 


925.8 


5.97 X 10-3 


3.94 X 103 


7.629 X 10-^^ 


Ne2+ 


36.02 


0-1 


399 


1.15 X 10-3 


7.20 X 10^ 


6.337 X 10-1^ 
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The critical electron densities for the transitions [ricr = Aui/kui, where kui is the elec- 
tronic coUisional de-excitation rate coefficient) are taken from Mendoza (1983) and are given 
in the next to last column of Table 1. We assume that electronic collisions are more important 
than H collisions in exciting these fine-structure transitions. 

We ignore the J = — 1 36.02 fim transition of Nelll, which is an order of magnitude 
weaker than the J = 1 — 2 15.55 fim transition, and the very weak — 2 magnetic quadrupole 
transition (not listed in Table 1). The optically-thin approximation is valid because the 
emission comes from a surface layer with limited vertical column density < 10^^ cm~^. 
Deeper regions (A^h > 10^^ cm~^) do not contribute because they are either insufficiently 
ionized or too cool for the transitions to be coUisionally excited. 

The flux in a transition m— >/ of an ion with charge n (the ion label n = 1, 2 now becomes 
a line label, n = 1 for the 12.81 and n = 2 for the 15.55 fim transition) from a disk annulus 
between R and R + dR is, 

1 

dF^i{R- n) = 2'KRdR / dN P„(n) {A^ihvui)n , (3-1) 



/o 



where Pu{n) is the normalized population of the upper level, A^i and hvui are the usual A- 
value and energy of the transition, and dN = dNfi is the differential of the vertical hydrogen 
column density. At high densities, the populations are given by the Boltzmann distribution, 
whereas at low densities they reduce to the probability for a coUisional excitation to occur 
within the radiative lifetime of the transition. In our calculations, the excited levels of the 
Ne fine-structure levels are populated near-thermally at i? = lAU, but substantially sub- 
thermally at larger radii. The departure from a thermal population is conveniently expressed 
in terms of the factors, 

n„{ul) 

Cui = l + 3-2 

which make it possible to use Boltzmann-like formulae for the populations in this problem. 
Thus the population of the upper level of the Nell 12.81 yum J = 1/2 — 3/2 transition is 
given by, 

^^^^^ = 2C3/2-i/2eii22.8/T + i' (3-3) 
and the population of the Nelll 15.55 /xm is given by 

^"^^^ " 1 + 5/3 Ci_2 e925-3/T + 1/3 Co_i ■ ^^'^^ 

Equation 13-41 is an approximation based on ignoring the weak quadrupolar transitions, both 
coUisional and radiative, that connect the top J = and the bottom J = 2 levels of the Ne III 
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fine-structure triplet. The factor Cui is unity in tlie fiigli-density limit (thermal population) 
and proportional to the electron density in the low-density (sub-thermal) limit where every 
exciting collision produces a fine structure photon. 

The key factors in Equation 13-11 for the differential flux are the integrals, 

iV„(i?;l) = J dN Puil) xi = J dNPuil)^, (3-5) 

N^{R; 2) = dNP^{2) x^ = x^, dN P^{2) ^^1±3^^ (3.6) 

making use of Eq. I2-I2j Nu{R; 1) and Nu{R; 2) are the column densities of excited Ne ions 
that emit the dominant fine-structure lines of Ne II and Ne III at radius R. Thus the radial 
distribution of the flux is given by, 

dFM = ^ {Auihvui)n 27r RNu{R] n) dR, (3-7) 
where the numerical values of A^ihuui are given in the last column of Table 1. 

In Figure 4 we plot Nu{R;n), which determines the emissivity per unit area of the 
12.81yum line of Ne"^ and the 15.55/im line of Ne^+ through Equation 13-71 Two extreme 
thermal models are shown, one for X-ray heating dominant {ah = 0.01) and the other for 
mechanical heating dominant {ah = 1.0). The emissivities are generally strongly peaked at 
very small radii where the upper levels are near-thermally excited. The effects of the central 
peaking is suppressed somewhat due to the factor R in Equation 13-71 for the flux. 

At larger radii, the excitation is sub-thermal, and the emissivities are determined by 
the second and third powers of the electron density, as suggested by Equation 12-111 The 
electron density profiles decrease in magnitude roughly as the inverse square of the radial 
distance. For large R, the emissivity is further reduced by the impact of the decrease of both 
He and T for large column densities. The emissivity of the 15.55/im line is smaller because 
of the reduced abundance of Ne^+, due in large part to its destruction by charge exchange 
with atomic hydrogen. In addition, the characteristic quantity {Auihuudn is smaller by a 
factor of 1.75 for this transition compared to the 12.81//m line. The 12.81/im line of Ne^ is 
more sensitive to the heating model because it tends to be emitted deeper in the atmosphere 
where mechanical heating can extend the warm temperature transition region. 

The unresolved line fluxes can be obtained by integrating the emissivities in Figure 4 
over radius, including the extra factor of R in Equation 13-71 The results are given in Table 2, 
assuming a nominal distance for the model T-Tauri star of = 140 pc and a Ne abundance 
of 10~^. Approximately half of the flux of the 12.81/im line of Ne+ arises from within 6 AU, 
with the remaining half produced over the radial range 6 — 30 AU. 
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Table 2: Neon Ion Line Fluxes (erg s ^ cm ^) 



ah 


Nc+ 12.81 




Ne2+ 15.55 


1.00 


1.25 X 10- 


14 


6.46 X 10"^^ 


0.01 


6.22 X 10" 


15 


5.27 X 10-16 
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The flux of the 12.81 //m line for the two extreme heating models differs by only a factor of 
two, so wc conclude that a typical T Tauri star generates a flux of ~ 10~^^erg cm~^s~^; the 
flux of the 15.55 fim line is about 20 times smaller, ~ 5 x 10~^^erg cm~^s~^. The predicted 
value of the flux of the 12.81 fim line suggests that it may be observable with only moderate 
spectral resolution. The dust continuum emission for the D'Alessio generic T- Tauri star (at 
140 pc) near 13;um is 1.21 x 10^^^ erg s^^ cm~^ /xm~^. If we consider the Spitzer Infrared 
Spectrometer (IRS), with a spectral resolution of 600, the continuum flux within a spectral 
resolution element is 2.58 x 10~^^ erg s~^ cm~^, about 25 times larger than the theoretically 
predicted line flux. Thus a high signal-to-noise observation with the IRS should lead to 
a detection of the 12.81 /xm hne. The Ne^"*" 15.55 fim line would be much more difficult 
to detect. These lines would become easier to detect in sources with weaker continua or 
with observations made at higher spectral resolution. Because H atom collisions have been 
ignored, the fluxes in Table 1 should be considered as lower limits. 

This discussion of the observability of the Ne flne-structure lines has been based on 
a neon abundance of 10~^, or equivalently on a solar abundance ratio of neon to oxygen, 
Ne/0 = 1/6. Several lines of evidence have suggested that the Ne abundance is closer to 
Xnc = 2.5 X lO^''. Chandra X-ray spectroscopy of nearby stars with a wide range of the ratio 
Lhoi/Lx give a nearly constant Ne/0 = 0.4 (Drake & Testa 2005). For the few T Tauri stars 
where the same measurements have been made, two are in harmony with this result (BP 
Tau and TWA 5 with Ne/0 0.5 ) and one has an even larger value (TW Hya with Ne/0 
= 0.87 ± 0.13), as discussed by Drake, Testa & Hartmann (2005). Measurements of Orion 
B stars also give a larger Ne abundance than 10~^, x^e = 1-29 x 10"^ and a larger neon to 
oxygen ratio, Ne/0 = 0.25 (Cunha, Hubeny & Lanz 2006). The solar Ne abundance plays 
a key role in the continuing controversy between the results of helio-seismology and solar 
models (Asplund et al. 2004, Antia & Basu 2005, Bahcall, Basu & Serenelli 2005, Schmelz 
et a/. 2005, Delahaye & Pinsonneault 2006). Despite these uncertainties, our results suggest 
that disk Ne abundances anywhere in the above range of values are sufficient to produce 
detectable Ne+ emission. 

The spatial distribution of the emissivity in Figures 4 has consequences for the rota- 
tional broadening of the Ne fine structure lines emitted by protoplanetary disks. Assuming 
Keplerian rotation, the column density Nu{R;n) can be converted into a distribution func- 
tion P{v) for the rotational speed {P{v)dv oc Nu{R)RdR) such that P{v)dv is the fraction 
of emission within the speed interval dv, 

P{v) (xv-^N^{R{vy,n) (3-8) 

where R = GM^v^'^. When combined suitably with the distribution of turbulent velocities, 
P{v) helps determine the shape of the lines. Figure 5 shows the distribution function P{v) 
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for the Nell 12.81 //m line, where the abscissa is w — v/v{l AU) > 0. We see that P{v) has 
two components. One is a narrow distribution centered at 0.25 AU), which corresponds 
to the cut-off of the spatial variation of Nu{R; n) in Figure 4 beyond 16 AU. Then there is a 
long tail extending to the maximum rotational speed, which arises from the peaking of the 
spatial distributions in Figure 4 at small R. Ignoring turbulent broadening, the rotational 
lineshape function for inclined disks will have a double peak, associated with the cut-off of 
the emissivity at large radii, plus extended line wings associated with the concentration of the 
emission at small radii. In principle, this analysis can be used to obtain information about 
the spatial distribution of the emission from measurements of the lineshape. Assuming that 
the turbulent velocities are of order lkms~^, as suggested by the measured CO hneshapes 
of T Tauri stars (Najita et a/. 2006), moderately-high spectral resolution observations with 
~ 3 X 10^ could confirm the expectations for the Ne lineshapes based on Figure 5. 



4. Discussion 

In the last section, we concluded that the Ne"*" 12.81/xm line should be detectable from 
a typical T Tauri star with a moderate resolution spectrometer like the Spitzer IRS. Indeed, 
we have reports from colleagues of the detection with this instrument of unusually strong 
12.81/im emission from a few T Tauri stars (llaria Pascucci 2006, Dan Watson 2006 private 
communications). These are probably the brightest Ne fine-structure emission line sources; 
the emission from most others is likely to be fainter. This line is also observable from the 
ground, and attempts to observe it with a large telescope at high spectral resolution would 
be of great interest. Not only would the line to continuum ratio be increased, but it might 
become possible to obtain information about the line shape, as discussed at the end of §3, 
that would help locate the source of the emission in the disk. 

We have argued that the detection of the Ne ion fine-structure fines would provide clear 
support for the role of stellar X-rays in determining the physical properties of the atmosphere 
of YSO disks, the region most accessible to observational study. The Ne lines complement 
the detection of the ro- vibrational transitions of CO and the UV fluorescence of H2 (reviewed 
by Najita et al. 2006). The latter transitions are associated with material at disk radii within 
a few AU and with vertical column densities Ah ^ lO^^cm^^. Our calculations show that 
the emission of the Ne fine structure lines arises from a larger range of radial distances, out 
to 10 — 15 AU, and that the lines are formed mainly at vertical column densities in the range 
Ah ~ 10^^ — lO^^cm"^. Thus the Ne fine-structure transitions are complementary to the 
lines of CO and H2 in probing both the warm temperature region above the cool molecular 
layer and a region that extends to relatively large radii that, in our own solar system, would 
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range from the location of the terrestrial planets to beyond the giant planets. By contrast, 
both the UV transitions of H2 and the CO ro- vibrational transitions are generally restricted 
to the inner few AU of disks (Najita et al. 2003; Najita et al. 2006), and the IR transitions 
of H2 are likely to probe a narrow range of radii (Pascucci et al. 2006). 

The measurement of Ne hnes in the disks of low-mass YSOs is also of interest because 
of the special chemistry of neon. From the disk modeling perspective, Nc has the advantage 
that the complexities of molecule formation and destruction do not enter. It is unlikely 
that significant amounts of Nc will freeze out on dust grains, although some of the neon 
trapped in meteorites has been interpreted in terms of adsorption (e.g.. Swindle 1988). 
Thus the detection of the Ne fine-structure lines, in conjunction with model calculations, 
might eventually provide information on the abundance of Ne in YSO disks, which serve as 
observable analogs of the primitive solar nebula. Because the Ne lines arise in a relatively 
thin layer at high altitudes, by themselves they provide only weak constraints on the total 
gas mass of protoplanetary disks. However, when taken together with other diagnostics, such 
as CO and H2, detailed modeling of several species might provide important information on 
disk gas mass. 

The calculations reported in this paper are limited to the fine-structure lines of just 
one element and one type of disk, the typical T Tauri disk as formulated by D'Alessio 
et al. (1999). In addition to extending the range of application of models such as ours (see 
also the work of Gorti & Hollenbach 2004), attention needs to be given to the underlying 
physics and chemistry of disk modeling, some of which is quite uncertain and affects the 
accuracy of model predictions. 

The emission of Ne fine-structure lines from relatively cool YSO disks arises from the 
ejection by hard X-rays of the K-shell electrons of moderately heavy ions, followed by the 
Auger effect. This process is enhanced as the mass of the target ion is increased, and it surely 
occurs for other abundant elements such as Si, S, and Fe. Our analysis of Ne fine-structure 
lines shows how important the recombination cascade from higher to lower ion stages is in 
determining the abundance of the first few (and most abundant) ions. Sulfur serves as a 
good example. Gorti and Hollenbach (2004) and Hollenbach et al. (2005) have suggested 
that the SI 25.2 and 56.3 ^m. fine-structure lines, among others, are potential diagnostics of 
intermediate-age YSO disks. Hard X-rays will produce ions as high as S^"*", and the highest 
ions are expected to be destroyed by a chain of fast charge transfers with atomic hydrogen, 
as has been confirmed by a high-level theoretical calculation by Stancil et al. (2001) of the 
reaction S''^ -|- H — S^+ + H+, which is fast down to low energies. However, the situation 
for 8^+ and 8^+ remains problematic. Early calculations of S^+ H — > S+ H+ (Butler 
& Dalgarno 1980, Christensen & Watson 1981) suggest that the rate coefficient is small at 
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low temperatures. Recent calculations of charge exchange (Bacchus-Montabonel 1998; 
Labuda et al. 2004) arc restricted to keV energies. The recent study of the reaction, H"*" 
+ S — >^ S"*" + H, by Zhao et al. (2005) provides important information for the ionization 
balance between S and S"*", but theoretical and experimental studies of and S^^ charge 
transfer with H at low energies (< 1 eV) are needed before the emission of the S 1 and S 111 
fine-structure lines can be calculated reliably for X-ray irradiated disks, e.g., along the lines 
of the analysis for neon in this paper. A similar conclusion applies to other heavy elements 
for which low-energy charge-transfer cross sections with H are lacking for the first few ions. 

It is important to recall the discussion of low-energy Nc ion charge exchange in §2. In 
particular, we have relied on 25 year old theoretical estimates in adopting a rate coefficient 
/c2 = 10~^^cm^s~^ for the charge transfer from Ne^+ to H. Because the detectabihty of 
the Nelll 15.55 //m line depends on this coefficient, it would be important to have more 
theoretical work on this reaction. We also have no information on the rate coefficients for 
the excitation of the Ne ion fine-structure levels by atomic hydrogen. Although few in 
number, existing calculations (e.g.. Launay & Roueff 1977 & Barinovs et al. 2005 for H 
-|- C^) lead to rates ~ lO^^cm'^s^^. If this order of magnitude is appropriate for the Ne 
ions, then H atom collisions would be important in that part of the disk atmosphere where 
< 0.01 and T is still large enough to excite levels with excitation temperatures of order 
1000 K. Our flux estimates may also be affected by the lack of understanding of other basic 
microscopic processes, not to mention the many astrophysical uncertainties such as disk 
structure, dynamics and thermodynamics. We have considered the possibility that the Ne 
ions are destroyed by molecules and PAHs. It is known from room temperature experiments 
(Anicich 1993) that Ne"*" docs not react strongly if at all with H2 and CO, molecules that 
arc predicted to occur at reduced levels in the warm upper atmosphere. Depending on 
their charge and abundance, PAHs might conceivably pay a role in destroying Ne ions and 
thus reduce the strength of the fine-structure lines. Negatively charged PAHs can neutralize 
Ne ions with large reaction rate coefficients (e.g., Tielens 2005). Both capture and charge 
transfer are possible reactions with neutral PAHs. Using rate coefficient estimates in the 
literature (Bohme 1992, Tielens 2005) and a maximal PAH abundance of 10~^, we find that 
these neutral reactions are less important than destruction by radiative recombination. Of 
course, there is also the strong possibility that the PAH abundance in YSO disk atmospheres 
is reduced by the strong stellar X-ray flux or by the ions that the X-rays generate. It would 
appear that the largest uncertainty in the underlying microphysics of our analysis of the 
neon ion abundances and their flne-structure line emission is in the magnitude of the weak 
but not insigniflcant charge transfer to atomic hydrogen. 

One astrophysical uncertainty that affects our estimates of the line fluxes is the GNI04 
choice of the X-ray luminosity and spectral temperature, Lboi — 2 x 10^°ergs~^ and kTx — 
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1 keV. The former is characteristic of solar-hke YSOs observed by Chandra in the Orion 
Nebula Cluster (Garmire et al. . 2000, Wolk et a/. 2005). This may well be an overestimate 
for the X-ray luminosity of lower-mass and older T Tauri stars. Based again on Chandra 
observations of the Orion Nebula Cluster (Preibisch et a/. 2005), the X-ray luminosity is 
measured to be smaller for lower-mass and older YSOs, and it is suppressed somewhat for 
accreting as opposed to non-accreting systems. Rather than use just a single temperature, it 
would be more appropriate to represent the X-ray spectrum with a two-temperature model 
of the type used to fit X-ray observations of well-observed stars. More generally, in future 
modeling of specific stars, the observed X-ray spectra should be used. 

We can use the approximations developed in §2, especially Equations 12-101 and 12-111 to 
identify some of the other astrophysical uncertainties in our flux estimates. The Ne levels 
are sub-thermally excited over most of the disk volume responsible for the emission of the 
fine-structure lines. Thus the local emissivity of the Nell 12.81 /im line is determined by n^, 
because the Ne"*" abundance usually tracks the electron fraction (Eq. 12-111) . Equation 12-101 
now tells us that the emissivity is determined roughly by 

j(Ne+)(xxNenHC(Ne)^^, (4-1) 

ai[T) 

i.e., by the product of the Ne abundance, the local density, the ionization rate (includ- 
ing attenuation), and a function of temperature that is essentially the Boltzmann factor, 
exp {—Tui/T). This last factor is one of several responsible for cutting off the emissivity at 
large column densities (and large radii). Equation 14-11 suggests that the Nell 12.81 /im line 
strength will be larger for younger and more active T Tauri stars, which have larger densities 
and ionization rates. However, this does not mean that the line will be easier to detect from 
such stars, because their dust continuum emission will also be larger. 

Flaring of the disk also plays a role in that the more flaring the smaller the attenuation 
of the X-rays. The GNI04 is flawed in this respect because it uses the D'Alessio et al. (1999) 
density distribution, whose scale height is determined by the dust temperature. In a thermal- 
chemical model that treats hydrostatic equilibrium self-consistently, we would expect an even 
more diffuse atmosphere than the one used in this paper. Furthermore, Table 2 shows that, 
without mechanical heating, the flux Nell 12.81 /im line is reduced by a factor of two. In 
other words, the emissivity is affected by the poorly understood processes that heat the disk 
atmosphere. Clearly, there are many possibilities for variation of the Ne fine-structure line 
fluxes, and some if not many T Tauri stars may generate fluxes smaller than those given in 
Table 2. 

In conclusion, we have modeled the physical properties of the disk atmosphere of a 
generic T Tauri star exposed to a strong stellar X-ray flux, and we have calculated the abun- 
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dance of Ne"*" and Ne^+ and the strength of their fine-structure hues. The estimated fluxes 
indicate that the Nell 12.81 /xm hne should be detectable with existing instrumentation. 
The observation of this line from warm disk atmospheres would provide strong evidence for 
the significant role of stellar X-rays on the physical and chemical properties of YSO disks 
and on the conditions for photoevaporation. 

The authors would like to thank Alex Dalgarno for advice on neon ion charge transfer 
reactions with atomic hydrogen, Phil Stancil for information about sulfur ion charge ex- 
change, Jack Lissauer for discussions on the potential role of neon in determining the gas 
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David HoUenbach for discussing the possible role of EUV radiation in generating Ne ions 
and especially for his many helpful comments on this manuscript. This research has been 
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Fig. 1. — Temperature profiles for a protoplanetary disk atmospliere based on GNI04. Tfie 
radial distance is 5AU, and the mass accretion rate is 10~^Mq yr~^. The abscissa is vertical 
column density in cm~^, and the ordinate is temperature in degrees K. The lower dashed 
line is the dust temperature for the D'Alessio et al. (1999) generic T Tauri star model. The 
upper solid curves are gas temperature labeled by the phenomenological mechanical surface 
heating parameter defined by Equation I2-I4t ah = 1 (mechanical heating dominant) and 
ah = 0.01 (X-ray heating dominant). The major chemical transitions occur between 10^^ and 
10^^ cm~^: significant levels of H2 and full formation of CO occur near ~ 1 — 2 x lO^^cm"^, 
and complete formation of H2 is achieved at Nn = 8.5 x 10^^ cm~^. 

Fig. 2. — Gas temperature profiles for the GNI04 disk atmosphere model for the case ah = 
0.01. The abscissa is vertical column density in cm~^, and the ordinate is temperature in 
degrees K. The general structure consists of a thin hot layer on top of a thick cool layer (near 
the midplane), with a warm transition layer in between. 

Fig. 3. — Electron fractions for the GNI04 X-ray irradiated disk atmosphere model for the 
case ah = 0.01. The abscissa is vertical column density in cm~^, and the ordinate is Xe, 
the electron abundance relative to total hydrogen. Deeper down than a vertical column 
A'h = 3 X 10^°cm~^, the electron fraction is < 10"^. 

Fig. 4. — The Ne fine-structure line emissivity per unit area, expressed in terms of column 
densities of the upper level from Equation l3-7[ calculated with the GNI04 disk atmosphere 
model for a neon abundance of XNe = 10~^. There is a pair of curves for each fine-structure 
transition, Nell 12.81/im and Nelll 15.55/xm. The upper curve is for ct/j = 1 (mechanical 
heating dominant), and the lower curve is for a = 0.01 (X-ray heating dominant). The 
emissivity is peaked at small R, and then drops off rapidly beyond R ~ 16 AU. 

Fig. 5. — The distribution of rotational speeds of the Nell 12.81 /im line emission. The 
abscissa is speed, normalized to the value at 1 AU (the smallest radius used in the present 
emissivity calculation), and the ordinate is the unnormalized probability. The peak is located 
near 0.25f(lAU), corresponding to the large- i? cut-off in Figure 4. The long tail at large 
velocities is associated with the strong emission at small radii. This distribution, convolved 
with the distribution of turbulent velocities, determines the line shape, which will be double- 
peaked and exhibit substantial line wings for inclined disks. 



-25- 




-26- 



10000 



ir~rT 



TTT 



TTT 



TT 



TTT 



I I I lllllj 



3000 — 



9 1000 
o 

II 

g 300 

•i-> 

u 
a 

6 100 — 

H 



30 — 




1 AU 



2 AU — 



_.5_AU 
.10.^U 

20 AU 
_30.^_ 

40AU 



10 
10 



I I 1 1 



i_LL 



J_LL 



J_LL 



J_J_L 



I I I 



18 



.19 



.20 



.21 



^22 



10 10 10 10 10 

Vertical Column Density 

f2.eps 



23 



10 



24 



-27- 




-28- 




1 2 5 10 20 50 

R (AU) 



f4.eps 



-29- 




f5.eps 



